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Abstract 
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these WIMPs might self-gravitate and eventually collapse forming a mini-black hole that eventually 
destroys the star. We impose constraints competitive to direct dark matter search experiments, 
for WIMPs with masses down to the TeV scale. 
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I. INTRODUCTION 



Observations of clusters of galaxies, rotations curves of individual galaxies, cosmic mi- 
crowave background anisotropies, and many other methods suggest the existence of dark 
matter. A possible realization of dark matter might be in the form of Weakly Interact- 
ing Massive Particles (WIMPs). A huge effort is being undertaken by experimentalists to 
directly detect WIMPs in underground or space experiments, as well as by theorists to 
incorporate them into viable theories beyond the Standard Model. 

The situation experimentally is still not clear, as the majority of the experiments have not 
detected WIMPs so far. Direct search experiments with Earth based detectors like CDMS [1] 
and Xenon [2] have imposed constraints on the WIMP-nuclei cross sections, assuming the 
local dark matter density around the Earth as inferred from the cosmological and other data 
(see e.g. Ref. [3] for the determination of the amount of dark matter from the WMAP data). 
On the other hand, DAMA experiment |1] claims dark matter detection with parameters 
that contradict other experiments if taken at face value. 

Given the still unclear picture regarding the nature of dark matter, it is of crucial im- 
portance to constrain as much as possible the WIMP candidates, including their mass and 
interactions. Several such candidates exist in the market depending on what theory be- 
yond the Standard Model one chooses, ranging from supersymmetry ^ E] and hidden 
sectors [3 [S], to Technicolor [S'-TIB]. The WIMPs can be classified according to their prop- 
erties, i.e. if they are produced thermally, if they are asymmetric [9l [161 E|) if they have 
spin- dependent or spin-independent cross section with the nuclei, if their collisions with the 
nuclei are elastic or inelastic [T8H22], and/or whether they are self-interacting [23H26] . 

Apart from direct searches, constraints on the properties of the WIMPs might arise from 
astrophysical observations as for example in [27]. Concentration of the WIMPs within stars 
can affect, under certain circumstances, the evolution of the latter, and/or products of WIMP 
annihilation within the stars could be directly or indirectly detected. The capture of WIMPs 
in the Sun and the Earth pSHSU] has been used to predict a possible signature for an indirect 
detection of dark matter based on neutrino production due to WIMP co-annihilation pH |32] . 
Constraints on the dark matter properties and the dark matter profile can also be imposed 
due to the effect of dark matter on the evolution of low mass stars [33l[3l] . and main sequence 
stars \35\, on possible gravitational collapse of neutron stars |36], and on the cooling process 
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of compact objects such as neutron stars and white dwarfs [371 - H2] . In particular, the authors 
of [30] have investigated under what conditions a neutron star can collapse gravitationally 
due to accretion of WIMPs, providing an upper bound for the WIMP masses given the 
local dark matter density and the time of accretion. Although the bound for the mass of a 
bosonic WIMP was low ~ 10 MeV, the upper bound for fermionic WIMPs was quite high 
~ 10^ TeV (and therefore not relevant for physics at the TeV scale). 

In this paper we investigate possible constraints that can arise from stars that accrete 
WIMPs during their lifetime and then collapse into a more compact object, white dwarf or a 
neutron star, inheriting the accumulated dark matter. Depending on the location of the star 
and the WIMP-nuclei cross section, it might be possible to impose constraints on the mass of 
the WIMP, excluding in some cases candidates that are lighter than TeV. Such constraints 
improve significantly the existing ones, and may become relevant for LHC physics. 

More specifically, we consider two different cases. In the first case, we examine the 
accretion of WIMPs with spin-dependent interactions with protons onto a Sun-like star. 
We deduce under what circumstances the accumulated WIMPs can trigger a gravitational 
collapse once the star has turned into a white dwarf. In the second case, a supermassive 
star accretes WIMPs which have a spin-dependent cross section with nucleons (protons 
and neutrons). A typical supermassive star of 15 solar masses lives about 10^ yr and then 
explodes forming a neutron star. Under certain assumptions, the WIMPs inherited by the 
neutron star from its progenitor will thermalize, sink to the center and, for some range of 
parameters, collapse further into a black hole. Thus, the mere existence of neutron stars 
might impose constraints on the mass or the cross section of the WIMPs. In all cases we 
will assume cross sections that are compatible with experimental constraints from the direct 
dark matter searches. 

II. ACCRETION OF WIMPS ONTO STARS 

The crucial parameter that determines the capture of WIMPs by a star is the WIMP- 
nucleon cross section. In this paper we will consider dark matter candidates that exhibit pre- 
dominantly a spin-dependent cross section with nucleons. Spin-independent cross section is 
highly constrained from direct dark matter search experiments, whereas the spin-dependent 
one is less constrained. The present constraint on the spin-independent cross section (nor- 



malized to a single proton) is roughly 3 x 10~'^^cm^ for a WIMP of a mass 100 GeV [T]. 
However, the rate of events is proportional to the local WIMP number density, and because 
for a fixed dark matter energy density the number density decreases inversely proportional 
to the WIMP mass, the present constraint can be written (for masses higher than TeV) as 



As one can see, the constraint becomes weaker at higher masses. 

The constraints on spin-dependent cross section are not so strict. The best upper limit 
on the spin-dependent cross section of WIMP with neutron is ~ 10^^^ cm^ |13], whereas 
the analysis for spin-dependent cross section between WIMP and proton gives a minimum 
upper limit by one order of magnitude higher |13] . The constraint for a spin-dependent cross 
section between WIMP-nucleon (proton or neutron) can be written as 



This constraint also weakens linearly in m at high masses. 

The difference between the spin-independent and spin-dependent constraints is due to 
several reasons. The first one is that WIMPs with the spin-independent cross section scatter 
coherently with the whole nucleus if their De Broglie wavelength is larger than the size of the 
nucleus. This condition is easily met on Earth based detectors. The coherence increases the 
cross section between WIMP and nucleus compared to the one of WIMP-nucleon roughly 
by a factor of A^^ where N is the number of nucleons composing the nucleus. For example, 
this gives an enhancement by a factor of ~ 73^ in the case of the Ge detectors of CDMS. 
In addition, form factors suppress further the spin-dependent cross section, making the 
resulting constraints even weaker. 

Let us now briefly summarize the accretion of WIMPs onto a star. Following [281 [371 E] . 
one can write the accretion rate as follows 



where pdm is the local dark matter density, v is the average WIMP velocity, M, and R are the 
mass and the radius of the star, Eq is the maximum energy of the WIMP per WIMP mass 
that can lead to a capture, and / denotes the probability for at least one WIMP-proton 
scattering to take place within the star. In this expression we have neglected relativistic 
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corrections (which are very small for regular stars) and possible motion of the star with 
respect to the dark matter halo. The latter can reduce slightly the accretion rate, but not 
more than an order of magnitude. Since we are not targeting particular stars like the Sun, 
we are going to present results for stars that do not move relatively to the dark matter 
halo. However, one can easily get the correct accretion rate in the case of a moving star by 
multiplying the accretion rate by (A/7r/2)erf(?7/?7), where r] = ^/3J2vq/v, with Vq being the 
velocity of the star. 

Let us first estimate Eq. The recoil energy produced by a WIMP-proton scattering, is 
within < T < Ampin/ {m + mp)^. Upon assuming that the distribution over the recoil 
energies is not very different from a uniform one, a typical scattering will produce a recoil 
of order ~ 2mp/m (for m » nip). In order for a WIMP to be captured (i.e. to become 
gravitationally bound) after one collision, it must lose at least the initial kinetic energy it 
had far out from the star. This leads to 

m n 

This is a conservative estimate because we have implicitly assumed that the collision will 
take place at the outskirts of the star (at the radius R), and not somewhere deep inside 
where the kinetic energy (and therefore the recoil energy) would be larger. 

For small cross sections, the probability / of at least one scattering of WIMP in a Sun-like 
star was estimated in Ref. |28] to be 

/ ~ 0.89— , (2) 

0"crit 



where 



cTcrit = -^-T- - 4 X lO^^^cm^ 



M \RqJ \Mq 

Mq and Rq being the mass and the radius of the Sun, respectively. Note that for another 
case of interest in what follows, a supermassive star of a mass M = 15Mq and radius 
R = 6.75i?0 the critical cross section is (Xcrit = 1-25 x 10~^^cm^. The probability / 
saturates to 1 for a > dcrit- 

We should also mention that in principle, the probability / may change as a function of 
time. As the star burns its hydrogen to helium, WIMPs interacting via spin-dependent cross 
section passing through the star meet fewer protons (hydrogen) scatterers to interact. He^ 
cannot interact through spin-dependent interactions with WIMPs. Therefore if we make the 



simple assumption that the star has converted all the hydrogen to helium at the end of its 
hydrogen stage, the probability / would drop to zero as WIMPs passing through the star do 
not find protons to scatter anymore off and be captured. It is natural to include this effect 
in the definition of the critical cross section. In this way one gets for the spin- dependent 

case 

where to is the star's lifetime and we have assumed that the initial composition of the star 
is 75% hydrogen and 25% helium, i.e. the one of Big Bang Nucleosynthesis. The critical 
cross section grows with time (the probability / decreases). 

The WIMP capture rate, Eq. ([T|, can be written in convenient units as follows, 

— (ofe) (^) (^) (I;) (i^) (- (^) 

where / is defined by Eqs. (|2|, and ([s]). Making use of Eq. ^ one can estimate the amount 
of accreted WIMPs during the lifetime of the star. A typical Sun-like star will burn hydrogen 
for 10 billion years before it becomes a red giant and later a white dwarf. A typical 15 Mq 
supermassive star burns first hydrogen to helium for 11.1 million years and then helium 
to carbon, carbon to oxygen etc. in a much smaller timescale (the last stage before the 
supernova explosion is the silicon burning lasting about 20 days). 



III. GRAVITATIONAL COLLAPSE OF THE WIMP-SPHERE 

Once the WIMPs are captured by the star, they start to thermalize through succes- 
sive collisions with the nuclei inside the star and after sufficient time are described by the 
Maxwell-Boltzmann distribution in the velocity and distance from the center of the star. 
The majority of WIMPs then is concentrated within the radius 

where T is the temperature of the star, pc is the core density and m is the mass of the WIMP. 
In the last equality we used typical values for a Sun-like star T = 1.5 x 10^ K, and pc = 
150 g/cm^. For a supermassive star of mass M = I5M0, the thermal radius is roughly an 
order of magnitude larger (using typical values T = 3.53 x 10^ K, and pc = 5.81 g/cm^ |15]). 



Depending on the mass and the cross section between WIMP and nucleon, WIMPs might 
or might not thermahze during the hfetime of the star. Since our constraints will depend on 
the WIMP thermalization, we estimate here the thermalization time. The thermalization of 
captured WIMPs can be divided into two stages: at the first stage the WIMPs oscillate in 
the star's gravitational potential crossing it twice per period. This lasts until the WIMP's 
orbit decreases to the size of the star. At the second stage, the WIMP moves completely 
inside the star on the orbit which shrinks to the thermal radius. 

Consider the first stage. Each time the WIMP crosses the star it has a chance to collide 
and lose some energy. The time between collisions At is given by half a period of WIMPs 
oscillation divided by the ratio of the WIMP cross section to the critical cross section. At 
each coUision the WIMP loses typically a fraction 2mp/m of its energy. Averaging over the 
WIMP trajectory inside the star (assuming for simplicity that the latter passes through the 
center) the typical energy loss is 



where M and R are the mass and the radius of the star, and r is the size of the WIMP's 
orbit. Dividing this energy change by At and expressing r in terms of energy gives the 
differential equation for the WIMP energy as a function of time. 



where — GMm/R is the binding energy of the WIMP at the star surface. Prom this 
equation the duration of the first stage is 



where eg = I-EqI/-^* is the ratio of the WIMP initial energy Eq to its binding energy at the 
star surface. Numerically eo ~ mp/m, and for a solar mass star and spin-dependent cross 
section we have 



The time becomes longer for larger masses and smaller cross sections. 

At the second stage, the average time between two successive collisions is At = l/{nav), 
where n is the number density of the nucleons, a is the WIMP-proton cross section, and 
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V = ^JlE jm is the average velocity of the WIMP. Therefore, the energy as a function of 
time is determined by the following equation 

= -navAE = -2V2pa - , (8) 
dt \m J 

where p is the matter density of the star. Solving this equation gives the time t2 needed to 

decrease the energy from E^ to Ef, 

~ ^pa \ ^fE~ )' 

Here the initial energy is of order ~ GM/R, while the final energy is determined by 
the final size of the WIMP cloud. In case of cooling to the thermal radius this energy is 
Ef = {3/2)kT, where k is the Boltzmann constant and T is the temperature in the core of 
the star. One can safely neglect the term depending on Ein in Eq. ^ since Ein ^ Ef. 

In case of cooling to some fixed radius (e.g. to the radius of a future white dwarf ~ 
4000 km, as will be of interest in what follows) the final energy is determined by this radius 
and the density of the star. Assuming typical parameters of a solar mass star, one gets 

t, = 0.15 yr ( ^) ( . (10) 

^ VTeV/ V 10-35 cm2/ ^ ^ 

Once inside either a white dwarf or a neutron star, WIMPs start to thermalize once 
again toward the much smaller thermal radius of a white dwarf (typically a few kilometers 
for a TeV WIMP) or of a neutron star (typically a few centimeters for a TeV WIMP). This 
process described by the same Eq. (|9]), but with different parameters. If a sufficiently large 
number of WIMPs has been accumulated, WIMPs may start self-gravitating and collapse 
gravitationally (in the absence obviously of a repulsive force between them). In this way 
the formation of a compact star (either a white dwarf or a neutron star) may trigger the 
collapse of the WIMP sphere into a black hole. 

The onset of the self-gravitating regime happens when the total mass of WIMPs inside 
the thermal radius becomes comparable to the total mass of the ordinary matter in the same 
region. This leads to the condition 

where we have substituted typical parameters of a white dwarf. For a neutron star the 
required number of dark matter particles is 5 orders of magnitude smaller. 
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The self-gravitating WIMP sphere may collapse into a black hole if the Fermi pressure of 
the WIMPs cannot counterbalance the gravitational attraction. The onset of the gravita- 
tional collapse occurs when the potential energy of a WIMP exceeds the Fermi momentum, 
and therefore Pauli blocking cannot prevent the collapse anymore. This happens when 



In the derivation of the above limit, we have considered that WIMPs are (semi)-relativistic, 
which is justified since once WIMPs self-gravitate themselves, they get closer and closer, 
building up a Fermi momentum that eventually corresponds to relativistic velocities. From 
the above equation we can deduce the number of WIMPs needed for the collapse to take 
place. 



where m is the WIMP mass and mpi is the Planck mass. 

IV. CONSTRAINTS ON DARK MATTER 

Having derived the accretion rate formula for a generic star and the amount of dark matter 
needed in order to form a mini-black hole, we can proceed to the constraints that arise from 
the requirement that such mini-black holes are not created inside newly-formed white dwarfs 
and neutron stars. We consider two different cases, i.e. constraints on spin-dependent cross 
sections from white dwarfs, and from neutron stars. 

A. Constraints on Spin-Dependent cross section from White Dwarfs 

Dark Matter WIMPs can have purely spin-independent, or spin-dependent interactions 
with nuclei, or even both types at the same time. Due to the coherence effect, the spin- 
independent interactions are usually stronger than the spin-dependent ones. However, there 
are cases where the spin-independent cross section is either suppressed or absent. Such 
cases arise naturally in models where dark matter candidates have an axial coupling to 
gauge bosons. One characteristic example is Majorana particles. Majorana fermions scatter 
off nuclei without the enhancement mentioned earlier because the amplitudes of scat- 
tering on different nucleons add up incoherently. Since most of the nucleons within the 
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nuclei come in pairs of opposite spin, the WIMP interacts effectively only with the un- 
paired nucleons. However, a Majorana particle is its own antiparticle, allowing therefore for 
WIMP co-annihilation. The co-annihilation invalidates our constraints because it destroys 
the WIMPs before they collapse into a black hole, unless the annihilation cross section is 
extremely small. However, in such a case the spin-dependent elastic cross section should be 
parametrically equally small. 

The constraints we present in this subsection are valid for dark matter candidates which 
predominantly have a spin-dependent cross section without being Majorana particles. We 
constrain models of asymmetric dark matter with WIMPs that have axial couplings to gauge 
bosons like Z . Such candidates have been identified and studied in For example, Dirac 
fermions with predominant axial coupling to the Z boson have suppressed spin-independent 
cross section and dominant spin-dependent one. In asymmetric dark matter models of this 
type, although the annihilation cross section is not suppressed, annihilations are rare due 
to the asymmetry. In such cases we can impose constraints competitive to the direct dark 
matter search experiments. 

In this subsection we look at potential candidates that have a spin-dependent cross section 
with protons which is larger than the spin-independent cross section. We consider both types 
of cross sections compatible with the experimental constraints. In such a case, WIMPs will 
be captured by a solar mass star primarily due to their spin-dependent interactions. When 
the star turns into a white dwarf, most of the WIMPs which are inside the white dwarf 
radius will be inherited by the white dwarf. 

In order to form a black hole, the accumulated WIMPs have to cool further. The spin- 
dependent interactions inside the white dwarf are suppressed because the latter is composed 
predominantly of spin-zero nuclei like He^, C^^ and O^^. However, since the white dwarfs 
are much denser than the ordinary stars, much smaller interactions are sufficient for the 
successful WIMP thermalization. As we will show below, these interactions may be provided 
by a small spin-independent component in the interaction and/or by a small admixture of 
nuclei with non-zero spin such as for example C^^ or O^^. 

Consider as an example the case where there are spin-independent interactions in addition 
to the spin-dependent ones. In this case the dark matter particles are accumulated in a solar 
mass star mostly due to their spin-dependent interactions, but once the star collapses to a 
white dwarf made of spin-zero nuclei, they thermalize due to collisions of spin-independent 
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FIG. 1. Constraints on the spin-dependent WIMP-proton cross section. The constraints follow 
from the existence of old white dwarfs in globular clusters with a core dark matter density of 
10^ GeV/cm'^ (straight solid line) and 10^ GeV/cm^ (straight dashed line). We also show the 
constraints from direct searches (solid curve). 

nature. We do not have to assume that WIMPs interact also with neutrons. It is sufficient 
to have spin-dependent interactions with protons only. In fact, the experimental constraints 
on spin-dependent cross section of WIMPs with protons is even less stringent, the strongest 
bounds coming from CDMS |13], PICASSO |17], and KIMS [48j. The constraint is approx- 
imately given by 



at WIMP masses above 1 TeV. Note that stricter constraints on the spin-dependent cross 
section between WIMP-proton from the ICECUBE collaboration |19] do not apply in our 
case since the constraints are based on annihilation of WIMPs, and we constraint candi- 
dates with spin-dependent cross section but of asymmetric nature, which effectively have no 
annihilations because the antiparticles are simply not present anymore. 
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In Fig. 1 we present the resulting constraints on the spin-dependent cross section, as a 
function of the WIMP mass. As one can see, these constraints are competitive with the 
direct ones at WIMPs masses of 1 TeV or less, depending on the dark matter density at the 
location of the observed white dwarf. The constraints cannot be extended to lower masses 
because no matter what the cross section is, not enough WIMPs can be accumulated to 
trigger the collapse. 

Let us now estimate the spin-independent cross section required for thermalization of the 
WIMPs inside the white dwarf. Substituting typical white dwarf parameter in Eq. ^ one 
finds ^ 

Given that cold and old white dwarfs (more than several billion years) have been observed 
in the inner parts of globular clusters [SUj, the thermalization time may safely be as long as 
1 Gyr. Assuming for simplicity that the white dwarf is made exclusively of C^^, this gives 

.&>4xl0-(— ) (— ) (^j cm^ (14) 

where to is the maximum time we allow for thermalization (e.g. 1 billion years) and a^j 
is the spin-independent WIMP-carbon cross section. The constraint can be rewritten as a 
constraint for the WIMP-proton cross section as 

.l,>8xlo-(^)-(^)-'^(A-)-^.^ (.) 

taking into account that the spin- independent cross section WIMP-carbon is related to that 
of WIMP-proton as ctsi/'^si — (/^c/yUH)^62, where fi are the reduced masses of the respective 
WIMP-nucleus and 6^ is the coherence enhancement due to the 6 protons of the carbon. 

As we have already mentioned, a spin-independent interaction of WIMPs with nucleons 
or a small fraction of isotopes of carbon or oxygen with nonzero spin can play the same 



role in the thermalization of WIMPs. By inspection of Eq. (13), it is obvious that one can 
always trade the cross section for the number density. This means that, assuming a TeV 
WIMP, if thermalization can be achieved for a^j ~ lO^^^cm^, it can equally be achieved 
by a spin-dependent cross section a^j^ = lO'^^cm^ where C now represents C^^ with a 
relative abundance ; (713 _ i ■ iQ^i^. Several red giants have been seen with a ratio 
(ji2 . (ji3 = 4 ; X predicted by the equilibrium processes. This is because C^^ from triple- 
alpha production can meet with hydrogen in outer shells of the red giant leading to C^^. 
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In practice, this means that even with a 10^ GeV WIMP an abundance of C^^ as low as 
Qi2 . (ji3 = 100 : 1 is sufficient for the thermalization. Note that the WIMP-C^^ cross 
section scales as ~ (^sof^c"^ / f^p times other nuclear spin factors, so the abundace of C^^ can 
be even 10^ smaller compared to C^^. We should also emphasize here that C^^ has an excess 
of a neutron, and therefore we have assumed that the WIMP couples equally to protons and 
neutrons. If the WIMP couples only to protons but not neutrons, another isotope has to be 
considered, such as O^^ or A^^^. 

For our constraints to be valid, we have to make sure that the black hole formed inside 
the white dwarf can eat/destroy the star within at most 1 billion years (smaller than the age 
of the older white dwarfs observed in globular clusters). First note that once the WIMPs 
form a black hole inside a white dwarf, regular nuclear matter starts falling in it. This 
process has been considered in Ref. [5Tj. It has been argued that for black holes of a size 
exceeding the atomic size the accretion proceeds in the Bondi regime which is characterized 
by a quasi-stationary matter flow into the black hole. The total rate of matter accretion 
can be expressed in terms of the matter density and sound speed far from the black hole. 
The accretion rate is proportional to the square of the black hole mass, so the change of the 
black hole mass with time is described by the equation 

where Mq is the initial black hole mass and is the characteristic time scale over which the 
star is destroyed, 

t - ' 



Here Cg and pc are the sound speed and the density in the core of the white dwarf, respectively. 
Numerically, the time is 

where we have estimated the value = 0.03c which is consistent with and the core 
density p = 10^ g/cm^. The heavier is the WIMP, the smaller is the amount of dark matter 
accumulated, and the smaller is the mass of the black hole created when the WIMPs collapse. 
For black holes made of collapsing WIMPs up to masses 10^ GeV or slightly lower, the time 
it takes for the full destruction of the star (i.e. the time it takes for the black hole to eat 
the whole star) is less than a billion years. 



The grey "no thermalization" area in Fig. 1, (see Eq. ([7]) and the discussion following) 
corresponds to cases where the captured WIMPs in the progenitor did not have enough time 
to settle within the radius of a white dwarf, meaning that although most of them will be 
gravitationally bound to the white dwarf when it is formed, their orbits do not necessarily 
intersect with it, and therefore there is the danger of not thermalizing and collapse gravita- 
tionally. Although in principle in such a case WIMPs might not collapse gravitationally, the 
situation is far from clear and we leave this for a future study. This is because WIMP- WIMP 
interactions can redistribute the angular momentum of the WIMPs and might eventually 
lead a large fraction of them to be captured fast by the white dwarf. In addition if the 
white dwarf is part of a binary system, there is again the possibility of allowing WIMPs 
to intersect with the white dwarfs. We should emphasize that to the left of the grey area, 
there is no ambiguity about the outcome, since WIMPs have enough time to concentrate in 
a radius smaller than that of the would-be white dwarf, and therefore once the white dwarf 
is formed, they are already inside. 

B. Constraints on Spin-Dependent cross section from Neutron Stars 

Constraints on spin-dependent cross section can be imposed also directly by accretion 
of WIMPs into neutron stars [36] , or by accretion of WIMPs in a progenitor supermassive 
star that later collapses to a neutron star. We consider supermassive stars that accrete 
dark matter during the hydrogen burning stage. At other stages of the star evolution 
the accretion is negligible since the star is mainly composed of spin- zero nuclei like He^, 
(ji2^ qw^ Q^f,^^ g^j^fj g^lgQ these stages are much shorter. After the last (silicon) stage, the 
supermassive star collapses to a neutron star. The supernova explosion cannot blow out 
the accumulated WIMPs because the spin-dependent interaction with silicon is zero. After 
the formation of the neutron star, the accumulated WIMPs are located mostly outside 
the neutron star, simply because the thermal radius of a supermassive star where WIMPs 
concentrate during the hydrogen-burning stage is several orders of magnitude larger than 
the radius of a neutron star. Upon assuming that WIMP- WIMP interactions or a binary 
system can redistribute the angular momentum and randomize the velocity of the WIMPs 
trapped by the supermassive progenitor, after the neutron star's formation, the WIMPs are 
captured by it very quickly [41]. Inside the neutron star the WIMPs continue to scatter on 
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neutrons until coming to a thermal equilibrium with the star. We should emphasize that 
we consider asymmetric WIMPs with spin-dependent cross section with both protons and 
neutrons here. 

There is an issue to address in the first place which is why to consider constraints which 
are due to the WIMP accumulation in the progenitor of a neutron star and not directly those 
which arise due to capture by the neutron star itself. The latter constraints have been first 
considered in [36] . For the spin- independent cross sections compatible with the experimental 
limits, a neutron star can accumulate a similar amount of WIMPs as a supermassive star. 
This is because the spin-independent cross sections consistent with the experiments are 
many orders of magnitude smaller than the critical cross section of a massive star, but not 
necessarily smaller compared to the critical cross section of a neutron star. Given the fact 
that massive stars live (and accrete dark matter) for a much smaller time compared to 
neutron stars, the latter are more efficient in accumulating WIMPs overall. However, in the 
case of spin-dependent interaction, experimental limits on the cross section are comparable 
to the critical cross section of a supermassive star. In such a case, a massive star can 
accumulate more WIMPs in the 10 million years of the hydrogen-burning stage than a 
neutron star would capture in 10 billion years. Therefore, the constraints resulting from the 
former process are stronger. 

Several pulsars have been detected in globular clusters close to (or inside) the core [5^ . In 
addition neutron stars older than 10 billion years have also been observed (as for example the 
B1620-26 in M4). In Fig. 2 we present the constraints on the spin- dependent WIMP-nucleon 
cross section that arise from the non-gravitational collapse of WIMPs inside neutron stars 
in dark-matter-rich environments such as the globular cluster M4. The grey rectangular 
area corresponds to exclusion due to direct neutron star accretion of WIMPs for roughly 
10 billion years (with pdm = 10'^ GeV/cm^) similarly to what previously has been derived 
in [361 [391 HI]. The area appears as a rectangular due to the fact that it excludes cross sections 
down to the critical cross section of a neutron star (~ 10~^^ cm^ [H])- The triangle area 
in purple is excluded (if one assumes the redistribution of WIMP velocities we have already 
discussed) due to WIMPs that have been accreted by the progenitor during its lifetime, 
and collapse in the core of the neutron star once they are "sucked" inside the neutron star. 
As one can see from the figure, the existence of a neutron star in the core of a globular 
cluster can impose in principle constraints competitive to those from the Earth based dark 




FIG. 2. Constraints on the spin-dependent WIMP-nucleon cross section that follow from the 
existence of neutron stars in globular clusters. The grey area is excluded by direct accretion of 
WIMPs (for 10 Gyr) in an old neutron star with local dark matter density 10^ GeV/cm'^. The 
purple area shows potential exclusion upon making an assumption regarding the distribution of 
WIMP velocities (see text), with the solid (dashed) line corresponding to local dark mater density 
10^ (10^) GeV/cm^. We also show the exclusion area by direct dark matter search experiments. 

matter search experiments for WIMP masses above 3 to 10 TeV depending on the local dark 
matter density. Our constraints become better than the direct ones when the mass increases 
further, since at larger WIMP masses fewer particles are needed for the collapse, whereas 
for the Earth based experiments it simply means fewer events in the detector. 

In order to make sure that our constraints are valid, we have to estimate the time it takes 
for the WIMPs to collapse once the neutron star is formed. As we mentioned earlier, the 
time it takes for the WIMPs accumulated by the progenitor to be captured by the neutron 
star is negligible. Once inside the neutron star, WIMPs have to lose energy via collisions in 
order to concentrate in the center and collapse. This time has been estimated in [4lJ and 
for all relevant cases is at most of the order of a year. In addition, we have checked that 
for the whole relevant range of WIMP masses, self-gravitation of the WIMPs sets on even 
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before they all concentrate within the thermal radius of a neutron star. This means that 
the collapse is even faster. 

We should emphasize that if a mechanism for randomizing the velocities of the WIMPs 
is operating (and therefore our constraints are valid), we exclude parameter space (purple 
triangle) that is not excluded by the Earth based experiments or by the gravitational collapse 
of WIMPs accumulated by the neutron star itself. Such a constraint can also be to some 
extent more robust than the one derived previously because no assumption is made regarding 
the age of the neutron star (as long as there is enough time for the black hole to destroy the 
star) . On the contrary, in the constraints that come from direct accretion of WIMPs in the 
neutron star, one has to know with enough precision the age of the star in order to estimate 
the amount of WIMPs accreted. 

V. CONCLUSIONS 

We derived constraints on the spin-dependent cross section of asymmetric fermionic dark 
matter WIMPs based on the existence of white dwarfs and neutron stars in globular clusters. 
Our constraints are competitive to direct dark matter search experiments, excluding a large 
parameter space of cross sections and masses as low as TeV (or slightly lower than TcV). 

In the case of white dwarfs, we were able to exclude a range of spin- dependent cross 
sections and WIMP masses because for these parameters, WIMPs that have been captured 
during the lifetime of the progenitor have enough time to concentrate within the core of the 
star that is inherited by the white dwarf, and eventually collapse gravitationally forming 
a black hole that destroys the star. This constraint is robust in the sense that it depends 
only on the local dark matter density of the globular cluster and no other hypothesis. We 
demonstrated that asymmetric WIMP candidates with only spin- dependent interactions 
with masses even lower than one TeV, trapped during the lifetime of the progenitor can 
easily thermalize inside the white dwarf due to expected small abundances of isotopes of 
carbon or oxygen that carry spin. 

If we now relax the strict condition of thermalization, i.e. to assume that WIMPs are 
gravitationally trapped by the progenitor but arc not necessarily confined within the radius 
of a white dwarf, our constraints can be extended to higher masses and lower cross sections. 
However, in this case we have to make an extra assumption of a mechanism that redistributes 
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the WIMP velocities in order for WIMPs which are on orbits around the white dwarf to 
intersect with it. Such a mechanism can be possibly provided by binaries or WIMP- WIMP 
interactions. Further investigation of this possibility is needed. 

In the case of neutron stars, we exclude an area of WIMP masses and cross sections 
due to direct accretion of WIMPs in the neutron star. The constraints depend only on the 
local dark matter density of the globular cluster and the age of the neutron star. Upon 
assuming an extra mechanism of WIMP velocity redistribution, extra parameter space may 
be excluded down to the TeV scale. 
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